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Abstract 

Starobinsky described an inflationary scenario in which quantum correc- 



tions to vacuum Einstein equations drive the inflation. The quantum cos- 
| mology of the model is studied by solving the Wheeler-DeWitt equation. A 

connection between uncertainty requirement, randomness in initial states and 
curvature fluctuation is studied with a Schrodinger-type equation through a 
time parameter prescription. The result obtained is applied to understand the 
5_i ■ decoherence mechanism in quantum gravity in the Starobinsky description. 

a i 
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I Introduction 

The Starobinsky scenario describes inflation through self-consistent solution of vac- 
uum Einstein equations, 

- = -8nG <T ixu >, (1) 

where < > is the expectation value of the energy momentum tensor. With a 
metric in Robertson- Walker form , 

ds 2 = dt 2 - a 2 {t)da k 2 , (2) 
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where da^ is the metric on a unit three sphere for k — 1, and the quantum correc- 
tions take a particularly simple form for free, massless, conformally invariant scalar 
fields: 

< T, u >= h + k 3 ^H, U) (3) 

where k\ and k 3 are constants and 

(!)^ = 2i W - 2^ V + 2i2i2^ - \g, v R 2 (4) 
2 1 1 2 

(3)^ = RffR^a — -RRfj, u — -g^iyR^ Rar + -^9nvR (5) 

The characteristic features of the model are as follows: 

(i) The model has a trace anomaly. Identifying 

SttG < T, u >= H- 2 + l -M - 2 (6) 

we get 

8ttG < T v v >= H - 2 (\r 2 - R VIJ R V °) - M~ 2 R.^ (7) 

though < T v v > vanishes for the classical conformally-invariant fields. 

(ii) The model shows de-Sitter solutions 

a(t) = Ho^coshiHot), k = +1 (8) 
a(t) = a exp(H t), k = (9) 
a(t) = Ho" 1 smh(H t) k = -I (10) 

These solutions describe an inflationary phase driven entirely by quantum correc- 
tions. 

(iii) Studying instability it has been found [1] that the solutions (8) - (10) are unstable 
under small perturbations and lead to a matter dominated universe at late times 
with a(t) oc £3. 

(iv) The curvature fluctuations are found to be 

5R 2 GH 2 

(-5") — ■ ( n ) 

Ro 7T 

The present work investigates the quantum cosmology of the Starobinsky model. 
Vilenkin [2] attempted this problem in order to study tunneling from "nothing" to 
the Starobinsky inflationary phase. The quantum analysis is then used to calculate 
the curvature fluctuations and assume a tunneling ansatz for the classical evolution 
of the model. 

In view of recent trends in solving the Wheeler-DeWitt equation [3-6], we take up 
study of the Starobinsky description again to understand the curvature fluctuations, 
quantum force driving the Starobinsky inflation, initial conditions for the classical 
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evolution from quantum phase. We also study the quantum to classical transition 
through the mechanism of decoherence. 

It is worthwhile pointing out that, at the semiclassical level, the Starobinsky 
inflationary scenario has been criticized [7] on the ground that the inflationary solu- 
tions are not perturbatively expandable in the parameter of "quantum corrections" 
terms in the equation of motion or Lagrangian. There is also some question of treat- 
ing the high order terms in the Lagrangian on an equal footing with the Einstein 
terms. Even then, studying the quantum cosmology of the Starobinsky description 
from the standpoint of quantum to classical transition would help us rethink the 
criticism of the model. Leaving aside this fact, it is instructive to look at the model 
as a toy example, to understand the current boundary condition proposals, decoher- 
ence mechanism and also the origin of quantum force in the early universe. Another 
motivation for studying the Starobinsky description is due to its inbuilt mechanism 
leading to spontaneous nucleation of the universe. Among the various types of in- 
flationary models, the choice of the inflaton potential has no microscopic origin, i.e. 
the choice of potential is not dictated from particle physics phenomenology (except 
Coleman- Weinberg potential). For this reason, the warm inflationary scenario [8] 
has been a crazy inclusion among the various types of inflationary models. In view 
of this fact, studying of the quantum cosmology of the Starobinsky model would 
be interesting, at least from the viewpoint of quantum to classical transition with a 
view to understand the initial conditions of the pre-inflationary era. 

In section II we obtain an approximate form of action for Starobinsky's descrip- 
tion thereby obtaining the Wheeler-DeWitt equation, and discuss some results of 
Vilenkin [2] for comparison with our work. In section III we discuss the time prescrip- 
tion in quantum gravity to reduce the Wheeler-DeWitt equation to a Schroedinger 
form containing time. The initial conditions are then treated in the light of our 
discussion. In section IV, we mainly concentrate on the decoherence mechanism to 
understand the quantum to classical transition and this constitutes the main text of 
the present paper. In this paper we work with wormhole-dominance proposal [9] as a 
boundary condition to the Wheeler-DeWitt equation and find that Gaussian ansatz 
is the corresponding initial condition for the Schrodinger- Wheeler-DeWitt equation. 
This correspondence is a breakthrough in the present trend of investigation along 
these lines [5,6]. Lastly, we end up with a discussion in section V. 



II Starobinsky Description 

The evolution equation for the scale factor is obtained from (1) ,(4), and (5) for 
k — and is given by 

TT 2 

H 2 {H 2 -H 2 ) = -^ 1 {2HH + 6H 2 H - H 2 ) (12) 
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In (12) H = | and the dot denotes the derivative with respect to cosmological time, 
and 

H = (8nk 3 G)-^, M = (48^0) (13) 

Equation (12) gives de Sitter solution H = H . Equation (12) for the condition 
H < H 2 , H « HH reduces to 

H 2 - H 2 = 6(ff,)H (14) 



L o 



with the solution 



iJ = iJ tanh(7--^-) (15) 

where 7 = | ln( ^" g ^~ g °^ ). Long inflation requires 

M 2 << 6H 2 . For H « H Q the term proportional to H A in (12) can be neglected. 

2HH + QH 2 H - H 2 + M 2 H 2 = 0, (16) 

with an approximate solution 

sin M n r 



rr 4 2 M t, 

H ~ — cos 

3t V 2 ' 



1 - 



Changing the origin of the time coordinates M Q t — > (M G t + a), we get 



a(t) = const, t 3 



i+ O smMot 



(17) 



(18) 



The expansion rate averaged over the oscillations period is 



ff = I < 19 > 

2 

so that a oc t » corresponds to the matter dominated universe. 

We have carried out [10] numerical integration of (1) with a contribution from 
particle production. We found that for a late time behaviour showing a radiation 
dominated (and also for a matter dominated) evolution, the initial state emerges 
mostly with a deSitter phase. Solutions (out of various solutions) with a(t) = a m i n 
at t — greatly effects both the initial and late time behaviour. This numerical result 
suggests that the classical universe sits for infinite time at the turning point a = a m i n , 
but quantum mechanically it enters into the disallowed region a < a min = H^ 1 and 
tunnels out again in the classical region with a curvature fluctuation ~ — ^ 
H 2 . To understand the quantum entanglement in the classical region, we need a 
quantum analysis in minisuperspace description. We need an action for the quantum 
description to get the Wheeler-DeWitt equation. For Starobinsky's description, (1) 
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or (12), no closed form of action is available. However for M Q << H 2 , a closed form 
of action is obtained in the form 

^^I + ^ + T^t' (20) 

with 

L(R) = ——(#+-% + — In— ) (21) 

V ; 16ttG V 6M 2 R Rj y ! 

where R Q = 12H 2 . The curvature fluctuation is 

( S_R )2 2 
{ R ) ~ ^|^ (jRo) | 



7T 



(22) 



as expected. 

To obtain the Wheeler-DeWitt equation through canonical quantization, one 
needs to determine the Hamiltonian H to obtain the equation 

= (23) 

In view of higher order terms like R 2 and R 2 In ■§- one gets 

S = J C(a,a,a)dt (24) 

in which the second derivatives a cannot be removed from (24). To eliminate a one 
writes 

S = I C(a,a,R,R)dt (25) 

and takes 

S = 2tt 2 / {L(R)a 3 -p[R- Qa' 2 (l + a + ad)] }dt . (26) 
Varying S with respect to R, we find 

P = a 3 L'(R) . (27) 

Putting (27) in (26) one can remove d through integration by parts. Introducing 
new variables q and x instead of a and R, 

q = Ho a(^ (28) 
^ = ^ln(^) (29) 
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and using the condition M 2 << H 2 one obtains the Wheeler-DeWitt equation with 
the replacement P q — > — i-§-, P x — > — iS- as 



<9 2 1 <9 2 



<9g 2 g 2 (9a; 2 



- V(q,x) 



ip(q, x) — 



In (30) 



2^ 



ffos) = A-V(1-(T + ^W) 
%(2x + e^-l) 



fix) 



where 



2H n 



67T 



We make the transformation Q = in (30) and get 

<9 2 1 <9 2 



dQ 2 Q 2 dx 2 



-Q 2 (l-Q 2 U(x)) 



ip = 



where 



= A(l -^(rr)) . 



(30) 

(31) 
(32) 
(33) 

(34) 

(35) 
(36) 



Vilenkin obtained a solution of (35) in the region q > q x , where q = q x ~ 1 + ^^(x) 
is the turning point, 



i> = **\Vo{qT k exp {-^(1 + f(x)) [l - t(q 2 - g 2 /i 2 - if 2 } + ^} 
where 

K(g) = A-V(l-g 2 ). 
We take g >> 1 in (37) and find 

V> ~ 7T^A5g _1 exp{-^-(l + fi 2 (x)) \l -iq 3 (l - fi 2 (x)f 2 ] + ^} 

3A ^ J 4 

1 



(37) 
(38) 



(ttA) V 1 exp {- — (1 + n 2 {x)) [l - *g 3 (l - fi 2 (x)f 2 } } . 



For x << 1, 



|?/>| 2 oc exp 



47TX 2 

GH ^ 



If one temporarily sets aside the normalization restriction on ■?/>, 



< 



5R X 2 o GH n 2 

>= 4 < x 2 >' 



2tt ' 



(39) 



(40) 



(41) 



in qualitative agreement with the classical result. The fluctuation in expansion rate 
is 

. 5H 15R ,GH 2 i 

6 ^ir = 2R-r { ^r )2 - (42) 

This 5 determines the duration [see (15)] of the inflationary phase in the model 

u~(^l- (43) 

Taking H Q ~ 0.7m p and using (42) and (43) we get 5 Q ~ 0.14, t* ~ 8jfa. Using 
8irki = 1.8 (restriction from SU(5) model) and the value of M Q from (13), t* ~ 
4 x 10 10 if This value of is more than sufficient to solve the horizon and flatness 
problem. The results quoted so far are standard and will be needed to understand 
much of what follows. 



Ill Reduction to Schrodinger Form 



Recent trends in quantum cosmology [3-6] suggest that the solution of (35) should 
be written as 

iP(Q, x) = exp [iS(Q)]$(Q, x) (44) 

in which WKB action S depends only on Q. However the philosophy behind this re- 
duction though differs, most of the workers obtained Schroedinger equation through 
a prescription on time variable. Apart from the interpretation, almost everybody 
remains silent about the normalization, and consequently also about the probabil- 
ity concept. These two concepts are of utmost importance for the interpretational 
framework of the wavefunctional ip(Q,x). We have shown [9] that if one takes 
contribution from wormholes in WKB ansatz, the wavefunctional becomes normal- 
izable. When boundary conditions corresponding to 'tunneling' and 'no boundary' 
proposals are introduced, one recovers the respective wavefunction. We have shown 
elsewhere that if one introduces the probabilistic concept through a continuity-like 
equation which also gives a time parameter prescription, the wavefunctional ip(Q,x) 
separates into the form (44) provided the classical Hamilton- Jacobi equation remains 
satisfied. We have thus a dynamical content to the form (44). We report briefly 
here the results only [11,12]. We put (44) in (35) and use 



and we find 



.<9$ 

'~dt 



d_ 

dt ~ ' 
1 d 2 



}_dS_d_ 
QdQdQ 



1 



2Q3 dx 2 + Q 2 U ( x ) + 2Q 



(45) 



(46) 



Remembering that in (35) remains multiplied by a factor Ti 2 andS in (44) ~ |r, 
we neglected terms which are of the order of h 2 . It may be noted that our 
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prescription does not violate unitarity which is obtained in the other prescription. 
The assumption (45) is consistent with (35) as well as with the Einstein classical 
equation. Equation (46) is the time dependent Schroedinger equation in quantum 
gravity. 

To understand the salient features of the Starobinsky description, let us proceed 
with the solution of (46) with Gaussian ansatz 

$(Q,x) = N(t)e-^ x \ (47) 
Substituting (47) in (46) we get coupled equations 

j - lnW= n +0 . A+ i ,48) 



With the anstaz 



q 3 m 



Q, = -iQ zV - (50) 

y 

S + 3§*-^» = 0. (51) 
Introducing conformal time coordinate rj according to dt = Qdrj, one finds 

y + 2 ~Qy — jp~ y = ^ ' 

where prime denotes a derivative with respect to 77. With the view that the classical 
model has a de Sitter solution , we take Q(rj) = where 77 runs from —00 to 0. 

Equation (52) now reads 

rj H Q r] 2 

which is solved by 

y = r? 3/2dz v /9/4+m 2 > (54) 

where we take m 2 = In the limit m 2 « | (which is usually assumed to be 

satisfied in the inflationary model and in our description) such that 



one finds from (49) 



^9/4 + m 2 ~ (3/2 + m 2 /3) (55) 
In the conformal time, the expression (50) reads 



>.y' 



n = -iQ 2 - (56) 

y 
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so that using (54) and (55) one finds 

n = -im 2 V\^-. (57) 
3 

As Q is imaginary the state (47) is not normalizable. To study decoherence one 
needs a real part in (57). We take up decoherence in the next section. 



IV Decoherence mechanism 

Decoherence is a mechanism through which we understand how the classical world 
arises from a quantum wavefunction of the universe. In decoherence, quantum 
interference effects are suppressed by the averaging out of microscopic variations 
not distinguished by the associated observables. In standard quantum theory this is 
referred to as the 'collapse of the wavefunction'. It is a formidable task to forbid the 
occurrence of linear superposition of states localized in far away spatial regions and 
induce an evolution agreeing with classical mechanics. In the context of quantum 
gravity the situation further complicates due to the absence of "time" because the 
Wheeler DeWitt equation when compared to Schrodinger equation 

ih^ = H1>, (58) 

gives a timeless character to the wavefunction of the universe. Though we recov- 
ered the form (58) with a reduced Hamiltonian given in (46), it must be ensured 
that during the quantum to classical transition , none of the successful quantitative 
predictions of the inflationary scenario for the present day universe is changed. An 
important aspect in this direction is to choose an initial condition for eqs. (58) or 
(35) such that basic input "inflation" remains undisturbed in the description. With 
respect to the Wheeler-DeWitt equation, the wavefunction if) is described with some 
boundary condition proposals namely (i) Hartle-Hawking [13] (ii) Vilenkin [14] and 
(iii) Wormhole-Dominance proposal [9] . The latter one is recently proposed by us. 

Now we will show that an initial adiabatic ground state with a Gaussian form is 
a suitable choice provided the quantum cosmological initial conditions correspond 
to the wormhole dominance proposal, or at least to the Hartle-Hawking proposal. 
This would also justify the correct choice of the boundary conditions. Using (57) in 
(48), we obtain for large Q 



N = N Q exp 



(59) 



The constant Nq will be evaluated through the wormhole dominance proposal. The 
wavefunction now reads, corresponding to (46) or (58), 



■ip = iV exp 



iQ 3 V\ , .M*V\Q 3 2 



(60) 
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We write the exponent in (60) as 

iQ 3 V\ 



S 



(1 - 

3 1 2Ef ] 



M 2 

Q 2 Ki - 



3/2 



A(l - ft* 2 ) 



(61) 



since jfix 2 « 1. Comparing with (33) we write (61) 



S 



3A(1 - n 2 (x)) 



3A(1 - n 2 (x)) 



[\Q 2 (1 



3/2 



3/2 



[g 2 (l -//(*)) 



--{l + ^x))[q 2 {l-^{x))-l 



3/2 



(62) 



The last step will be justified from the Wheeler-DeWitt equation itself since q >> 1. 
The wormhole dominance proposal considers incorporating the repeated reflections 



from the turning points q = o and q = q. 



r to contribute to N a as 



Here 

Evaluating (63) we find 



Hence 



cxp 



(l-Ai 2 Gr))S 

expS(q x ,0) 
~ 1-exp [2S(q x ,0)\ ' 

S(q x ,0) = \S\^. 

= exp{^(l + // 2 (rr))} 
° (l-exp{2/3A(l + /i 2 (*))}) 

£(1 + ^(x)){l-i(q 2 (l-^ 2 (x))- 1)1} 



(l-exp{2/3A(l + At 2 (x))}) 
We continue (66) in the region q 2 {l — yU 2 ) < 1 to get 

e^{Ul + ^ 2 {x))\l-{l-q 2 + q 2 ^ 2 f 2 ]} 



(63) 
(64) 

(65) 

(66) 



(l-exp2/3A(l + /i 2 (x))) 



(67) 



If we leave aside the denominator of (67), the wavefunction (67) corresponds to the 
Hartle-Hawking wavefunction. 

Thus the initial condition for the Schrodinger- Wheeler-DeWitt equation (46) 
turns out to be that at an early time (near the onset of inflation), the modes are in 
their adiabatic ground state due to wormhole dominance and this serves as the initial 
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conditions on the Wheeler-DeWitt wavefunction. Unlike the claim by some authors, 
it also establishes that Hartle-Hawking initial conditions also provides inflation and 
serve as a seed for decoherence to suppress the interference for quantum to classical 
transition. One important feature of (67) is that at q — > 0, (67) reduces to 

* 6XP{ ^ } (68) 

V (l-exp{2/3A(l + /i 2 (:r))}) ' 1 ] 

characteristic of the Hartle-Hawking wavefunction. As is evident from (34), A is a 
small quantity, and hence, considering the denominator of (68), we find 

|^| 2 ^ e -£(i+f 2 (*)) 

This gives the curvature fluctuation as in (41). 

The steps leading to the form (68) from (60) indicates that at the very early stage 
(within the turning points) the quantum superposition principle has been effective 
such that the time has lost its meaning. The interference between e lS and e~ tS is 
maintained through the wormholes that act as a driving quantum force. We have 
shown elsewhere [15] that without any reference to the Wheeler-DeWitt equation one 
can also recover eq. (46) exactly using only the classical Einstein-Hamilton- Jacobi 
equation with a directional time derivative 

* = A (69) 

dt QdQdQ' y ! 

i.e., the scale factor (gravitational field) itself acts as time. Just at the classical 
turning point, time begins to flow and the effect of the quantum turning point (at 
which the wormhole contributes) only survives through N(t) in the Gaussian ansatz 
i.e., the states e tS and e~ tS begin to decohere. It has been guessed by some authors 
[16] that some sort of boundary conditions at small scales may lead to quantum 
effects in the vicinity of the turning point. The 'wormhole dominance' proposal 
exactly establishes this aspect, keeping the necessary coinage for inflation. If we are 
to bridge the solution (47) with (67), the reduction 

S = S (Q)+S 1 (Q,x) (70) 

with 

■0 = exp [iSi(Q,x)\ 

is an unavoidable fact. Kiefer [5] named this reduction as the 'relevant' and the 
'environmental' degrees of freedom. The solution (67) refers to the ground state 
wavefunction. The higher degrees of freedom ('higher multipoles') when taken into 
account, eq. (53) modifies to 

y"--y' + (n 2 -^- 2 )y = o. (71) 
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Evaluating Q along the lines of [5] we find 

n 2 2 M 2 3 

and the decoherence factor responsible to suppress interference terms is written as 

exp(-£> ± ) = exp(-^Tr^-) 

= exp£n 2 A 2 . (73) 

n "i? 

We choose the adiabatic vacuum state with positive frequency as the solution 
of (71) or (53) and this acts as initial vacuum state. This initial state evolves at 
7] — > +oo according to a WKB form, 

y = * e ir«WW + i-i/VW , (74) 



where a and (5 are the usual Bogolubov co-efficients. For the WKB state (74) one 
finds 

n ~ lu + — (75) 
2w 



l~2 ^ ■ m 2 QQ' 

= ^n 2 - m 2 Q 2 - %—— — — , 76 

v 2{n 2 — m 2 Q 2 ) 



with m 2 = ^-A. Hence 



n 2 

* 8 V(n 2 -m 2 g 2 ) 3 ' 1 J 

Using 

1 00 2<f 

c _7rcot7rC = ?^e (78) 

and after some straightforward manipulation we find 



nm 3 QQ 2 
128 



7r . 9 ^, cot7rmQ 



1 + cot TxmQ) + — — —f- - 27r 2 cot nmQ(l + cot 7rm<5) 



mQ " (mQY 

(79) 

In the limit mQ » 1, one finds 

exp( _ fl±) . exp( _^ Hexp( _^, (SO) 

This result is also obtained by Kiefer [6]. For m ~ 100 GeV, H ~ 55km/secMpc, 
and Q — -f^" 1 , this is 

exp(-10 +43 ) ~ 0. 

Hence decoherence is efficient for large Q even in the Starobinsky description. This 
argument leads us to reconsider the Starobinsky description despite the criticisms 
labelled against it. 
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V Discussion 



The Starobinsky model provides a viable description both in the classically allowed 
and in the forbidden region. It is a model where one loop quantum corrections spon- 
taneously allow the universe to make the transition to a deSitter phase and then to a 
classical Friedmann stage with a small perturbation. It is a model of spontaneous nu- 
cleation rather than inflation-driven inflation. In our previous work [9] we achieved 
the advantage of normalized wavefunction with probabilistic interpretation but there 
still remains no certainty about the adoption of a particular boundary condition 
proposal. Within the perspective of the Schroedinger-Wheeler-DeWitt equation de- 
scribing the time evolution of the universe, we found that Hawking type boundary 
conditions definitely establish inflation (since we take Q = —^j^), as the character- 
istic feature of Starobinsky model, unlike the claim by some workers in favour of the 
tunneling proposal. The wormhole dominance proposal serves as an initial condi- 
tion for decoherence to effect quantum to classical transition, because the Gaussian 
ansatz for SWD equation (valid in the classically allowed region) leads, after con- 
tinuation, to the wavefunction according to the wormhole dominance proposal. We 
observe that the Starobinsky model (basically reminding us of a i? 2 -cosmology) is 
equivalent to an Einstein gravity plus a scalar field and this inflaton-like scalar field 
arises from one loop quantum corrections. It may be mentioned that quantum cos- 
mology with Hartle-Hawking boundary conditions for a model Einstein gravity plus 
a dialaton was studied by Okada and Yoshimura [17]. They briefly comment that the 
nucleation of the classical universe would be exponentially suppressed if one takes 
the tunneling proposal. The quantum nature of the early universe thus guarantees 
the universal validity of superposition principle. This in turn incorporates multiple 
reflections between the turning points. The wormhole dominance proposal takes this 
fact justifiably defining a normalization constant on the basis of superposition prin- 
ciple favouring probabilistic interpretation rather than the concept of 'conditional 
probability'. In our previous work [9] we have shown that the normalization constant 
prescription in the wormhole dominance proposal is equivalent to the contribution 
of wormholes leading to an interpretational framework for wormhole picture (lack of 
which created a confusion about its introduction), apart from Klebanov and Cole- 
man's arguments, as a driving quantum force in the early universe. The ensemble 
of quantum universe models would thus be a realistic situation in quantum gravity. 
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